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ABSTRACT 
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This work deals with an instability arising from the shear-induced migration of particles in 
dense suspensions coupled with a dependence of viscosity on particle concentration. The analysis 
summarized here treats the inertialess ( Re = 0) linear stability of homogeneous simple shear flows 
for a Stokesian suspension model of the type proposed by Leighton and Acrivos (1987). Depending 
on the importance of shear-induced migration relative to concentration-driven diffusion, this model 
admits short-wave instability arising from wave- vector stretching by the base flow and evolving into 
particle-depleted shear bands. Moreover, this instability in the time-dependent problem corresponds 
to loss of ellipticity in the associated static problem {Re = 0, Pe = 0). 

While the isotropic version of the Leighton-Acrivos model is found to be stable with their 
experimentally determined parameters for simple shear, it is known that the stable model does not 
give a good quantitative description of particle clustering in the core of pipe flow (Nott and Brady 
1994). This leads to the conjecture that an appropriate variant on the above model could explain 
such clustering as a two-phase bifurcation in the base flow. 

INTRODUCTION 

The general theme of the present work is material instability, whose recent history, from the 
1950s onwards, mainly resides in solid mechanics. There one encounters a wide variety of ’’strain- 
softening” phenomena giving rise to spatial localization of strain in the form of shear bands or 
” necks”. By comparison, there is a much smaller literature dealing with related instabilities in 
fluid-like materials, such as non-monotone viscous stress in particulate suspensions and polymeric 
fluids or dissipative clustering in rapid granular flows. 

A recent report (Goddard 1996) contains a brief review and fairly extensive bibliography 
whose common thread is the ” short- wavelength” instability of homogeneous deformations and the 
associated ” change of type” in the underlying field equations, e.g. loss of static ellipticity or of 
dynamic hyperbolicity (Joseph 1990). 

The present paper deals with a class of material instabilities in which diffusional scalar 
transport plays a crucial role. Prior related studies involve the effects of heat conduction on 
thermally-softening media, where mechanical dissipation (having no direct counterpart here) gives 
rise to phenomena such as bifurcation in liquid flows or ’’adiabatic shear bands in solids. 

MODEL AND ANALYSIS 

We consider a neutrally-buoyant ’’Stokesian” suspension (c/. Leighton and Acrivos, 1987) 
with stress: 

T = 2ty(0)D-pl (1) 
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and particle flux: 


(2) 


j = “ 7) +!/(</>, 7 ) V 7 } 

where <p denotes particle volume fraction, 

D := I{Vv + (Vv) r }, 7 := |2tr(D 2 )} 5 

and v is mixture velocity. The quasi-static mechanics (Re = 0 ,Pe> 0) is governed by the above 
constitutive equations, the balances 

V • T = Vp with V • v = 0, ( 4 ) 

and 

d t (j> + v • V<£ + V • j = 0, ' ( 5 ) 


together with suitable boundary conditions (compatibility with the base state and/or regularity of 
perturbed states in unbounded domains). 

For planar infinitesimal perturbations of the uniform base state: 

<4 0) = 7 (0) y and {^ 0 ) , v<°>, V 7 (0) , V<£ (0) } = 0, ( 6 ) 

the standard stream-function cum Fourier representation in the x-y plane ^(k-r k t) 
$(k x ,ky,t), must satisfy V V ' h 

A: 4 ^ + a(kl - &J)$ = 0 

and 

d t $ - k x dk v <$ = -k 2 {v(kl ~ k l)V + «$} 

where we let 

and adopt the change of notation 


(7) 

( 8 ) 

(9) 


7 (0) t 
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’ 7(0) 


^(0) 

^( 0 ) 


V 


( 10 ) 


The stream function ^ is readily eliminated from (7)-(8), which are easily seen to be elliptic for 
k > ua and otherwise hyperbolic, with the possibility of singular surfaces. 


Wave- Vector Shearing 

With 5 ' eliminated, (7)- ( 8 ) become 

d t $ - k x d kv <b = ct { 1 c }$ ( 11 ) 

where 

o{k} = -k*{K -*£(%- ktf} (12) 

Were it not for the term k x d ky , representing a wave-vector shearing (Thomson 1887, Savage 1992, 
Wang et al. 1996) that becomes dominant for large time, the preceding relation would immediately 
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provide the dispersion relation for mode k which, incidentally, is identical with that given by the 
simplified 1-D treatment in the concluding section of Nott and Brady (1994). 

We note that wave- vector stretching in an arbitrary isochoric homogeneous flow, with Vv = 
jT =cons t. and tr{L} = V • v = 0, can be concisely represented by the respective transformations 
between material (’’embedded”) coordinates (x,k) in physical space and ’’reciprocal” space: 


x = F *x, with (d t )x = ( dt) x + (Lx) • <? x 

(13) 

and „ x „ 

k = F r k, with (dt)^ = (dt)k — (L k) • 

(14) 

where 

F(t) = e Lt 

is the physcial-space deformation gradient. 

For the simple shear at hand, the transformation (14) gives k(k,t) as 

(15) 

k x — k x and ky k y k x t 

(16) 

and converts (11) into the o.d.e. 

dl °f = a{k(k,t)} 
at 

(17) 

where <r{k} is given by (12). . 

While (17) can be integrated subject to (16), asymptotic stability is determined by its 

limiting form for large t, and there are two cases to consider, k x ^ 0 and k x = 0, for which, 

respectively, one easily finds from (16)- (17) that 

k 2 t 3 kyt 

log$ ~ (va k) ~(i/o — k) — 

(18) 

and 

log$ = {va - K)k y t 

(19) 


f or i _* oo. Thus, k < va implies short-wavelength instability in the form of a (” Kelvin- mode ) 
shear-band structure, with k y » k x , representing particle-depleted strata lying perpendicular to 
the y-direction. 


A Special Case 


To describe the particle distribution in a fully-developed rectilinear shear flow, Leighton 
and Acrivos (1987) employ (1) together with the 1-D form for particle flux : 



jy = -^ 2 { K ^ + K2a ^} 

(20) 

where 

O 

A1 

c- 

-2 

"C 

II 

Q 

(21) 


and where the coefficients K, > 0 and K 2 > 0 are constants, or only weakly dependent on <f> (and 
proportional to particle diameter squared). 
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Taking (2) as the obvious extension of (20) to 3D (i.e. ignoring shear-induced anisotropy), 
one can easily deduce from the above stability analysis that linear stability (against planar distur- 
bances) requires that 

k 2 > Ki, ( 22 ) 

plainly satisfied by the empirically assigned values of Leighton and Acrivos (1987), for which K 2 « 
2K\ over a relatively broad range of 0. 

Recently Nott and Brady (1994) have made the interesting observation that the stable form 
of (20) does not give a particularly good quantitative description of the the dense ’’plug- flow” clus- 
tering of particles in the core of pipe flow. While they propose a more elaborate constitutive theory, 
this author wonders whether a mere adjustment of the coefficients K U K 2 , to impart instability 
near the state of maximum packing <t > max , might not suffice to explain the above clustering as a 
flow bifurcation with two-phase, core-annular structure . 

CONCLUSIONS 

Depending on its coefficients, the above particle-migration model can exhibit instability in 
simple shear, the dominant linear mode being particle-depleted shear bands. While the empirical 
form proposed by Leighton and Acrivos is linearly stable, an unstable version might allow for a 
better description of particle clustering in pipe flow. 

It is plausible that the migrational instability considered here would be enhanced by non- 
Newtonian rheology, as in shear-thinning suspensions of colloidal or deformable particles, thereby 
accentuating the tendency towards clustering or phase separation. 

For the simple constitutive model of the present work, loss of ellipticity in the static problem 
is directly connected with time-dependent instability. This broaches an interesting and unresolved 
mathematical question as to the precise relation between material stability and static ellipticity for 
materals with long-range memory. 

Several of the above issues are the subject of further work in progress. 
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